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In  Phase  two. all  currently  available  data  involving  inherent 
uncertainties  in  some  material  and  geometrical  parameters  was 
assembled  in  a  probabilistic  framework  and  subsequently  compared 
to  the  probabilistic  predictions  of  PROFRANC  qualitatively  and 
quantitatively.  These  comparisons  were  shown  to  be  very  suc¬ 
cessful.  The  quantitative  comparison  was  performed  by  hypothesis 
testing,  which  is  a  mathematical  rule  deciding  whether  to  except 
or  reject  PROFRANC  predictions  using  experimental  data,  f  |  j 

Results  to  date  show  in  particular  that: 


1.  PROFRANC  can  predict  deterministically  load  to  fracture 
initiation,  initial  direction  of  crack  extension,  com¬ 
plete  trajectory  to  failure,  and  various  load- 
displacement  relationships  in  arbitrary  two-dimensional 
structures . 

2.  PROFRANC  can  accurately  predict  probabilistic  character¬ 
istics  of  fracture  phenomena  for  systems  including  uncer¬ 
tainties  in  geometrical  and  material  parameters.  This 
predictive  capability  includes  not  only  means  and  vari¬ 
ances  of  all  relevant  fracture-based  responses,  but  also 
their  full  probabilistic  descriptions.  Two  independent 
hypothesis  tests  namely,  Kolmogorov-Smirnov  and  Wilcoxon 
tests,  validated  these  descriptions. 

3.  It  is  possible  to  create  a  controlled,  scientifically 
designed  database  for  both  deterministic  and  probabilis¬ 
tic  response  of  two-dimensional  mixed-mode  fracture 
propagation.  Such  a  database  currently  cannot  be  found 
in  the  open  literature  and  has  been  critical  in  the 
validation  of  PROFRANC. 

4.  A  previously  unanticipated  source  of  uncertainty  was 
uncovered  during  the  experiments  and  simulations: 
variability  in  supposed  boundary  conditions.  This  source 
of  uncertainty  has  been  shown  to  affect  significantly 
trajectory  of  curvilinear  fracture,  and  can  have  a  pro¬ 
found  effect  on  the  reliability  of  fail-safe  design  of 
structures  sensitive  to  such  cracking. 
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1.0  INTRODUCTION 

This  document  is  the  final  report  for  the  research  project 
entitled,  ’’Probabilistic  Fracture  Mechanics".  Its  principal 
objective  is  to  demonstrate  validation,  via  recent  unique  exper¬ 
iments,  of  previously  reported  theoretical,  analytical,  and 
numerical  developments  on  this  project.  These  developments  are 
focused  on  the  creation  of  a  new  capability  to  predict  probabil¬ 
istically  fracture  initiation  and  complete  history  of  propagation 
under  general  mixed-mode  conditions  in  arbitrary  two-dimensional 
structures.  This  new  capability  integrates  classical  mixed-mode 
fracture  mechanics  and  advanced  concepts  of  reliability  within  a 
state-of-the-art  computational  mechanics  code,  PRO FRANC,  housed 
in  an  engineering  workstation. 

The  validation  efforts  reported  herein  concentrated  on  pre¬ 
dictions  of  load  to  fracture  initiation,  initial  angle  of  propa¬ 
gation,  total  trajectory  to  structural  failure,  and  load-versus- 
displacement  and  load-versus-CMOD  relationships.  These  predic¬ 
tions  had  to  be  compared  to  observations  which,  because  they  are 
absent  from  the  literature,  had  to  be  generated  within  this  pro¬ 
ject.  Consequently,  an  experimental  program  was  designed  to 
develop  a  broad  database  for  mixed-mode  fracture  involving  uncer¬ 
tainties  in  material  and  geometrical  parameters. 

Section  2  of  this  report  reviews  the  scarce  available 
experimental  results,  while  Section  3  presents  details  of  the 
experimental  program.  A  two-phase  task  was  undertaken  to 
validate  PROFRANC.  Phase  one  consisted  in  predicting 
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deterministically  the  outcome  of  a  subset  of  a  larger  number  of 
experiments  in  which  variability  in  geometry  and  material  param¬ 
eters  was  purposely  minimized.  This  subset  included  the  experi¬ 
ments  with  variable  initial  crack  length  and  crack  location.  The 
purpose  of  this  phase  was  to  verify  that  PROFRANC  could  predict 
nearly  deterministic  events  accurately.  Section  4  shows  that 
this  phase  was  highly  successful.  Moreover,  a  more  realistic 
structure  has  been  tested  and  deterministic  validation  performed. 
The  structure  has  holes  of  variable  location  relative  to  initial 
crack  location.  Section  5  highlights  these  results.  In  Phase 
two,  all  currently  available  data  involving  inherent  uncertain¬ 
ties  in  some  material  and  geometrical  parameters  was  assembled  in 
a  probabilistic  framework  and  subsequently  compared  to  the  proba¬ 
bilistic  predictions  of  PROFRANC.  A  quantitative  comparison 
between  the  available  data  and  PROFRANC  preditions  was  performed 
by  hypothesis  testing.  These  comparisons  are  shown  to  be  suc¬ 
cessful  in  Section  6.  Furthermore,  this  section  includes  an 
immediate  practical  application  of  the  newly  developed 
technology. 
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2.0  PRELIMINARY  VALIDATION  EFFORTS 

2.1  Validation  of  Data  from  Literature. 

An  extensive  literature  search  has  been  conducted  on  experi¬ 
mental  data  for  crack  stability,  initial  direction  of  crack 
extension  and  crack  trajectory.  From  a  very  limited  number  of 
such  available  data,  a  reliable  set  of  nine  specimens  (Figure 
2.1)  was  found  [1]  for  a  preliminary  validation  of  our  theo¬ 
retical  model.  There  are  two  groups  of  specimens  in  this  set  of 
nine  corresponding  to  two  different  loading  conditions.  The 
groups  are  designated  as  800  series  consisting  of  six  specimens 
and  900  series  containing  the  remaining  three. 

Data  reported  in  this  study  show  that  geometrical 
parameters  such  as  initial  crack  length,  specimen  thickness, 
loading  angle,  and  material  fracture  toughness  KIC  exhibit 
statistical  variation.  Therefore  these  parameters  were  modeled 
as  random  variables. 

The  PRobabilistic  FRacture  ANalysis  Code  (PROFRANC) ,  devel¬ 
oped  at  Cornell  by  this  research  team,  was  applied  to  predict  the 
probability  of  load  at  fracture  initiation  and  of  initial  direc¬ 
tion  of  crack  extension.  Probability  distribution  of  the  frac¬ 
ture  initiation  load  was  in  excellent  agreement  with  the  exper¬ 
imental  results  (Figures  2.2  and  2.3).  The  experimentally 
measured  values  of  these  loads  are  very  close  to  the  peak  of  the 
corresponding  probability  density  curve  predicted  by  PROFRANC. 
Sensitivity  studies  show  that  fracture  toughness  is  the  critical 
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Load  lo  Ini  Halt  Crack.  PqlKNl 


Figure  2.3  Probability  density  curve  of  the  load  to  initiate  fracture  for  the 
900  series  specimens.  Experimental  results  are  also  shown. 
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parameter  for  fracture  initiation  load.  Measured  crack  initia¬ 
tion  directions  varied  by  about  two  degrees  from  predicted  direc¬ 
tions  (Tables  2.1  and  2.2). 

The  distribution  of  a  random  variable,  e.g.,  angle  of  or 
load  to  fracture  initiation,  predicted  by  PROFRANC  is  based  upon 
a  mathematical  model  derived  in  [3].  As  in  case  of  any  theoreti¬ 
cal  model,  this  model  also  relies  on  certain  assumptions  and 
approximations.  It  is  thus  essential  to  ascertain  qualitatively 
and  quantitatively  the  performance  of  the  model  by  comparing  the 
results  predicted  by  the  model  with  those  of  the  field  experi¬ 
ments.  In  this  section  we  have  studied  the  qualitative  compari¬ 
son  between  the  results  predicted  by  PROFRANC  and  laboratory 
experiments  from  literature.  For  completeness  of  the  model 
validation  study,  we  will  perform,  in  the  next  section,  quanti¬ 
tative  comparison  between  the  theoretical  and  experimental 
results,  referred  to  as  hypothesis  testing.  The  central  idea  of 
hypothesis  testing  is  as  follows.  We  have  a  random  variable, 
e.g.,  load  to  fracture  initiation.  From  the  theoretically  pre¬ 
dicted  (i.e.,  hypothesized)  distribution  of  this  load  we  can 
generate  its  samples.  On  the  other  hand,  we  can  generate  samples 
of  the  load  by  laboratory  experiments.  We  want  to  compare  these 
two  sets  of  samples  in  order  to  ascertain  whether  the  theoretical 
samples  belong  to  the  same  random  variable  as  that  of  the  experi¬ 
ment.  There  are  standard  methods  of  comparison  based  upon  sound 
mathematical  judgment.  We  shall  use  some  of  these  methods  as 
tools  of  our  hypothesis  testing. 


8 


Case 

800 

Series 

90% 

98% 

r.v.;  a,  ICjg 

(-67.54,  -67.56) 

(-67.54,  -67.57) 

A)  oty  Kit 

(-67.20,  -67.92) 

(-67.17,  -67.96) 

Table  2.1  Confidence  intervals  in  degrees  on  initial  direction  of  crack 

extension  for  the  specimen  in  Figure  2.1  using  the  parameters 
of  the  800  series,  Table  4,  [1] 


I 

I 

I 


Case 

900 

Series 

90% 

98% 

r.r.:  a,  Kje 

(-65.49,  -65.56) 

(-65.48,  -65.57) 

(kj  ty  CLy  K  fc 

(-65.10,  -65.91) 

(-65.04,  -65.96) 

Table  2.2  Confidence  intervals  in  degrees  on  initial  direction  of  crack 

extension  for  the  specimen  in  Figure  2.1  using  the  parameters 
of  the  900  series,  Table  4,  [1] 
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2.2  Hypothesis  Testing 

The  probability  density  curves  for  load  to  fracture  initia¬ 
tion  are  quantitatively  validated  in  this  section  by  two  estab¬ 
lished  statistical  tests  of  hypothesis.  Here  hypothesis  is 
referred  to  a  statement  about  the  probability  distribution  of  a 
random  variable,  namely,  the  load  to  fracture  initiation.  The 
test  of  hypothesis  is  a  rule  deciding  quantitatively  whether  to 
accept  or  reject  the  hypothesis  using  the  sample  data  of  the 
random  variable. 

2.2.1  Kolmogorov-Smirnov  Test  [33] 

This  procedure  uses  individual  sample  points  to  perform  the 
test  and  hence  is  effective  for  both  small  and  large  number  of 
sample  data.  The  test  compares  the  sample  (empirical)  cumulative 
density  function  (CDF)  of  the  random  variable  with  the 
corresponding  predicted  density.  The  test  consists  of  the 
following  steps. 

1.  order  n  samples  of  experimental  loads 
pql  -  pq2  -  *  *  *  -  pqn 

Sample  CDF:  FExp(pqj)  -  j/n,  j  =  1,  . . .  n. 

2.  Null  hypothesis,  HQ:  r.v  Pq  has  predicted  CDF  Fpred(pg) 
given  by  PROFRANC 

Alternative  hypothesis,  H1:  r.v  Pq  does  not  have  CDF 
FPred<Pq> 

3.  select  7,  the  level  of  significance  of  the  test 

4.  specify  the  rejection  region: 

R  ■  { D_.,v|  D__v  >  D  7  ) 

1  mc*.x  i  max  v  ' 
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Dn7  are  provided  in  statistics  Tables. 

5.  calculate 

Dmax  _  (lFPred(Pqi!  *  FExp<Pqi>l< 

l<i<n  ^ 

lFPred(pqi}  “  FExp(pqi-l^  I  * 

6.  if  Dmax  >  Dn7/  reject  HQ,  accept  H^. 

The  results  of  Kolmogorov-Smirnov  test  for  load  to  fracture 
initiation  are  given  in  Tables  (2.3)  and  (2.4). 

2.2.2  Wilcoxon's  Test  [34] 

Suppose  we  have  randomly  and  independently  drawn  n  observa¬ 
tions  from  the  population  of  experimental  loads  to  fracture  and  m 
observations  (m  >  n)  from  the  population  of  the  load  predicted  by 
PROFRANC.  Wilcoxon  test  checks  the  null  hypothesis,  HQ,  that  the 
two  populations  are  identical  against  the  alternative  hypothesis, 
H^,  that  they  are  nonidentical.  The  test  consists  of  the  follow¬ 
ing  steps. 

1.  Draw  m  samples  of  P  from  its  predicted  distributions 

H 

FPred^pq) *  Let  n  be  the  number  of  its  experimental 

observations . 

2.  Order  these  m  +  n  observations. 

3.  HQ:  r.v  Pq  has  predicted  CDF  Fpred(pq) 

Hx:  r.v  Pq  does  not  have  CDF  Fpred(pq) 

4.  select  7,  the  level  of  significance  of  the  test 

5.  specify  the  rejection  region: 

R  -  <w|w  < 

Wm  n7  are  provided  in  statistics  Tables. 
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Confidence  Limit,  7 

800  series) 

Test  Results 

i?3(900  series) 

Test  Results 

.01 

.618 

Accepted 

.828 

Accepted 

.05 

.521 

Accepted 

.708 

Accepted 

.1 

.470 

Accepted 

.642 

Accepted 

.15 

.436 

Accepted 

.597 

Accepted 

.2 

.410 

Accepted 

.565 

Accepted 

A n«=  i 


.22  800  series 

.548  900  series 


Table  2.3  Results  of  analyses  by  Kolmogorov-Smirnov  Test.  Random 
Variable:  Kjq 
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Confidence  Limit,  7 

DJ(  800  series) 

Test  Results 

£3(900  series) 

Test  Results 

.01 

.618 

Accepted 

.828 

Accepted 

.05 

.521 

Accepted 

.708 

Accepted 

.1 

.470 

Accepted 

.642 

Accepted 

.15 

.436 

Accepted 

.597 

Accepted 

.2 

.410 

Accepted 

.565 

Rejected 

.25  800  series 

.581  900  series 


Table  2.4  Results  of  analyses  by  Kolmogorov-Smirnov  Test.  Random 
Variable:  Kic,a,a,t 
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6 .  calculate 

W  =  X  Ri  ®  sum  of  ranks  of  experimental  observations 
among  m  +  n  samples. 

7.  if  W  <  Wm  n7,  reject  HQ,  accept  H^. 

The  results  of  Wilcoxon's  test  for  load  to  fracture  initiation 
and  initial  direction  of  crack  extension  are  given  in  Tables 
(2.5)  and  (2.6)  respectively. 

The  tables  show  that  the  predicted  distribution  of  load  to 
fracture  initiation  is  accepted  by  two  independent  tests  of 
hypothesis  for  almost  all  levels  of  confidence.  This  validates 
the  model  quantitatively  and  completes  the  preliminary  validation 
studies . 

Successful  validation  of  a  probabilistic  model  depends  upon 
the  number  of  experimental  data  of  the  random  variable  used  for 
the  test.  The  credibility  of  a  hypothesis  test  increases  with 
the  number  of  experimental  data.  In  the  above  studies,  the 
number  of  experimental  data  used  from  literature  was  limited. 

Our  extensive  search  shows  that  very  limited  amount  of 
experimental  data  are  available  in  literature.  This  necessitated 
the  additional  experiments  be  conducted  to  validate  current 
analytical  developments. 
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Table  2.5  Results  of  analyses  by  Wilcoxon  Test.  Random  Variable:  K|q 
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3.0  MIXED-MODE  FRACTURE  EXPERIMENTS 

3.1  Objective  of  Experiments 

The  very  limited  number  of  data  in  the  literature  necessi¬ 
tates  that  further  experiments  be  conducted  for  validation  of  the 
present  analytical  developments.  The  parameters  affecting  mixed¬ 
mode  fracture  initiation  and  trajectory  are  length,  inclination, 
and  location  of  initial  crack,  and  ratio  and  orientation  of  prin¬ 
cipal  material  toughnesses.  The  objective  of  the  experiments  is 
to  investigate  sensitivities  of  load  to  fracture  initiation  and 
crack  trajectory  due  to  uncertainty  in  the  selected  parameters. 
The  experiments  will  then  create  a  broad  base  of  experimental 
data  for  further  validation  of  our  analytical  developments. 

3.2  Background 

During  the  past  three  decades  many  experiments  have  been 
performed  to  study  the  problem  of  mixed-mode  fracture.  Erdogan 
and  Sih  [2]  performed  experiments  on  plexiglas  plates.  Each 
plate  had  an  inclined  center  crack  and  was  subject  to  uniform 
tension.  Williams  and  Ewing  [4],  and  Finnie  and  Weiss  [5]  tested 
the  same  specimen  configuration  with  different  materials.  Ewing, 
Swedlow  and  Williams  [6]  tested  plexiglas  panels  with  an  inclined 
edge  crack.  They  applied  two  types  of  loading:  uniform  tension, 
and  pure  bending.  Ingraffea  [7]  tested  plexiglass  and  rock 
plates  with  inclined  cracks  loaded  in  compression.  Banks-Sills 
[1]  tested  the  plexiglass  specimen  shown  in  Figure  2.1. 

Ingraffea  [8],  Arrea  and  Ingraffea  [9],  Iosipescu  [10],  Walrath 
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and  Adams  [11],  Sullivan,  Kao  and  Van  Oene  [12],  Barnes,  Kumosa 
and  Hull  [13],  Wang,  Goetz  and  Corten  [14],  Bocca,  Carpintteri, 
and  Valenti  [15],  Ballatore,  and  Carpinteri  [16],  and  Bazant  and 
Pfeiffer  [17]  tested  the  so-called  Iosipescu  shear  specimen  with 
different  types  of  materials  including  rock,  concrete,  composites 
and  polymers. 

In  most  of  the  reported  experiments,  only  the  initial  frac¬ 
ture  initiation  angle  and/or  failure  load  were  measured. 
Furthermore,  most  of  the  data,  such  as  initial  angle  of  crack 
growth,  was  reported  as  an  average  value.  This  does  not  allow 
one  to  make  any  statistical  analysis  of  previously  reported  test 
data.  Because  of  this,  and  because  of  interest  in  the  complete 
crack  trajectory  and  the  need  for  sufficient  variability  in  geom¬ 
etry  and  material,  many  additional  experiments  had  to  be 
performed. 

In  the  forthcoming  sections,  the  specimen  material  and  geom¬ 
etry  selection  is  described,  followed  by  a  description  of  speci¬ 
men  preparation.  The  experimental  setup  and  testing  procedure 
are  then  described.  This  is  followed  by  a  section  presenting 
some  typical  results. 

3.3  Selection  of  Specimen  Material  and  Geometry 

3.3.1  Selection  of  Isotropic  Material.  The  isotropic 
material  selected  for  use  is  commercial  plexiglas  (cyro  acrylic 
FF  plexiglas  MC)  manufactured  by  Rohms  Company.  This  material 
has  been  chosen  because: 
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(i)  It  is  a  good  approximation  to  a  homogeneous,  isotropic 
material . 

(ii)  On  the  scale  of  the  experiments  to  be  performed,  it 
behaves  in  a  linear  elastic  manner  and  obeys  princi¬ 
ples  of  linear  elastic  fracture  mechanics. 

(iii)  The  material  is  unaffected  by  changes  in  humidity  and 
temperature.  This  reduces  considerably  variability  in 
our  experimental  results  caused  by  unwanted  environ¬ 
mental  effects. 

(iv)  It  is  easy  to  machine  and  to  form  natural  cracks. 

(v)  It  is  transparent  so  that  it  is  easy  to  track  crack 
growth,  and  to  detect  accurately  fracture  initiation. 

(vi)  It  is  affordable,  so  many  specimens  can  be  tested. 

3.3.2  Selection  and  Preparation  of  Anisotropic  Material.  To 
study  the  effects  of  ratio  and  orientation  of  principal  material 
toughnesses,  it  was  necessary  to  find  an  affordable  material, 
isotropic  in  stiffness,  but  anisotropic  in  toughness.  After  an 
extensive  investigation  of  material  science  and  polymer 
engineering  literature  [18,19,20],  and  consulting  with  several 
researchers  at  Cornell,  it  has  been  found  that  by  hot-stretching 
plexiglas  above  its  glass  transition  temperature  (about  110 
degree  centigrade) ,  the  material  can  be  made  an: sotropic  in 
toughness  while  keeping  the  same  isotropic  elasticity  properties. 
Furthermore,  by  hot-stretching  the  material  to  different  strains 
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(from  0%  to  150%)  a  material  with  varying  toughness  ratios  (from 
1  to  About  2)  can  be  obtains i. 

To  create  the  specimens  for  this  project  the  plexiglas  was 
hot-stretched  at  a  temperature  of  120°C  to  an  elongation  of  about 
50%.  The  material  was  then  cooled  slowly  at  constant  length. 

The  cooling  rate  was  approximately  2°C/minute. 

3.3.3  Selection  of  Specimen  Geometry.  Figure  3.1  shows  the 
type  of  specimen  selected  for  testing.  It  consists  of  a  simply 
supported,  centrally  loaded  beam.  The  initial  crack  is  placed 
away  from  the  center  of  the  beam  to  produce,  under  the  given 
boundary  conditions,  mixed-mode  crack  initiation  and  curvilinear 
trajectory.  This  specimen  geometry  has  been  chosen  to: 

(i)  ensure  plane  strain  conditions  and  a  small  process 
zone  ahead  of  the  crack  tip  compared  to  crack  length 
and  other  characteristic  dimensions  of  the  specimen, 
based  on  principles  of  linear  elastic  fracture 
mechanics. 

(ii)  prevent  out-of-plane  buckling  and  minimize  maximum 
load. 

(iii)  guarantee  sufficiently  large  area  for  tne  crack 
trajectory,  and  to  enable  variable  geometry 
parameters. 

(iv)  ensure  that  the  specimen  may  be  easily  prepared  and 
loaded. 

A  similar  specimen  configuration  was  used  by  Jenq  and  Shah 
[21]  to  study  mixed-mode  fracture  in  concrete.  Note  that  the 
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selected  specimen  is  quite  similar  to  the  Iosipescu  shear  spec¬ 
imen.  The  reasons  for  not  using  the  Iosipescu  specimen  are: 

(i)  It  contains  two  initial  cracks  symmetric  to  each 
other.  In  practice,  one  will  always  initiate  before 
the  other.  This  will  make  the  specimen  response 
asymmetric,  and  difficult  to  analyze.  This  was 
pointed  out  by  Ingraffea  and  Panthaki  [22]. 

(ii)  the  selected  specimen  is  much  easier  to  prepare  and  to 
load. 

(iii)  the  selected  specimen  still  allows  mixed-mode  fracture 
if  the  initial  crack  is  placed  away  from  the  center  of 
the  beam. 

3.4  Preparation  of  Specimen 

Specimens  have  been  made  from  plexiglas  sheets  with  mechani¬ 
cal  properties,  provided  by  the  manufacturer,  shown  in  Table  3.1. 
However  the  key  properties  were  determined  directly  from  the 
material  used  in  the  tests. 

After  cutting  the  specimens,  all  edges  were  carefully 
machined  to  produce  dimensional  accuracy  of  +/-  1/32  inch.  At 
the  position  of  the  starter  crack,  a  notch  was  carefully  made 
using  an  automatic  bandsaw.  The  notch  length  is  accurate  to  +/- 
1/64  inch  and  orientation  was  accurate  to  +/-  1/4  degrees.  At 
the  end  of  the  starter  notch,  a  natural  crack  was  formed  by 
forcing  a  razor  blade  into  the  end  of  the  notch. 
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Property 

ASTM 

Method* 

Units 

Plexiglas 

MC 

Mechanical 

Tensile  Strength  ('/«"  specimen- 

0.2~/Minute) 

D638 

Maximum 

psi 

10.000 

Elongation.  Maximum 

% 

4.5 

Modulus  of  Elasticity 

PSI 

450.000 

Flexural  Strength  (Span  depth  ratio  16. 

O.r /Minute) 

D790 

Maximum 

psi 

15.000 

Deflection.  Maximum 

inches 

0.6 

Modulus  of  Elasticity 

psi 

450.000 

Compressive  Strength  (0.05‘VMinute) 

D695 

Maximum 

psi 

16.000 

Modulus  of  Elasticity 

PSI 

430.000 

impact  Strength 

0256 

Charpy  Unnotched 

ft.  lbs./ 

7.0 

’A"  x  1 "  section 

Rockwell  Hardness 

D785 

— 

M-90* 

Resistance  to  Stress 

Critical  Crazing  Stress 

ARTC  Mod. 

isopropyl  Alcohol 

Of 

psi 

•  1.300 

Toluene 

MIL-P-6997 

psi 

1.200 

Table  3.1  Mechanical  properties  of  plexiglas  provided  by  Rhoms 
manufacturing  company 
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Two  small  plexiglas  pieces  were  then  glued  to  the  specimen 
at  the  two  sides  of  the  notch  mouth.  A  clip  gage,  which  was  used 
to  measure  the  crack-mouth-opening  displacement  (CMOD) ,  was  then 
attached  to  these  pieces.  This  will  be  explained  in  more  detail 
in  the  next  two  sections.  Finally,  a  transparent  graph  paper  was 
glued  to  one  side  of  the  specimen  to  record  visually  and 
approximately  the  location  of  the  propagating  crack  tip  during 
the  test. 

3 . 5  Experimental  Setup 

The  experimental  setup,  shown  schematically  in  Figure  3.2, 
consists  of  a  new  MTS  load  frame,  a  control  microconsole,  a  data 
acquisition  system  with  A/D  and  D/A  converters,  a  microcomputer, 
and  a  microscope. 

An  MTS  control  microconsole  is  used  for  initial  control  of 
this  closed  loop  system.  It  is  used  prior  to  testing  to  zero  all 
the  channels  and  then  to  transfer  control  to  the  computer  program 
which  will  control  testing.  The  data  acquisition  system  contains 
both  A/D  and  D/A  converters  plugged  directly  into  the 
microcomputer.  An  interactive,  menu-driven  program  was  written 
to  control  the  experiments. 

Three  channels  are  monitored  during  the  experiments:  load, 
displacement  under  the  load,  and  crack-mouth-opening- 
displacement.  The  experiments  were  conducted  under  CMOD  control 
to  ensure  slow  crack  growth.  If  tested  under  load  control,  crack 
growth  would  always  be  unstable  and  dynamic.  On  the  other  hand, 


MTS  220  Kip  Testing  Frame 


Figure  3.2  Schematic  description  of  experimental  setup 
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if  testing  was  performed  under  displacement  control,  then  crack 
growth  might  or  might  not  be  stable  and  slow.  This  depends  on 
length,  inclination  and  location  of  the  initial  crack. 

Clausing  [23]  confirms  this  result  theoretically  based  on 
principles  of  linear  elastic  fracture  mechanics.  Figure  3.3 
shows  a  simply  supported  beam  with  a  center  crack  subject  to  two 
point  loads.  The  figure  reveals  that  there  is  a  transition  point 
from  stability  to  instability  when  one  decreases  the  initial 
crack  length.  The  same  conclusions  were  reached  by  Mai  and 
Atkins  [24,25].  Since  their  theoretical  results  were  only  valid 
for  Mode  I  tests  under  displacement  control,  the  applicability  of 
this  concept  to  the  mixed-mode  case  was  tested  by  performing  a 
series  of  stability  studies  using  the  finite  element  method.  It 
was  observed,  as  predicted,  that  length,  inclination  and  location 
of  initial  crack  are  the  key  parameters  that  control  specimen 
stability.  Under  CMOD  control,  the  selected  configuration  is 
unconditionally  stable. 

3 . 6  Test  Procedure 

Readings  of  load,  displacement  under  the  load,  and  CMOD  were 
performed  continuously.  The  CMOD  was  monitored  with  an  MTS 
extensometer  placed  to  measure  the  relative  opening  of  the  crack 
mouth.  The  extensometer  was  calibrated  over  the  expected  range 
of  CMOD  to  a  maximum  error  of  about  0.25%.  The  load  was  monitored 
through  a  load-cell,  mounted  directly  in  the  activation,  with  a 
capacity  of  5  kips  and  calibrated  to  a  maximum  error  of  less  than 
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0.1%  over  this  range.  Displacement  under  the  load  was  measured 
with  an  LVDT  mounted  within  the  load  actuator.  Although  this 
LVDT  has  been  calibrated  to  a  maximum  error  of  0.25%  over  its 
0.50  in  range,  its  reading  of  beam  displacement  was  relative 
since  support  settlements  and  load-line  displacement  is  also 
included  in  its  readings. 

The  specimen  was  loaded  monotonically  in  CMOD  control  using 
very  small  increments.  Prior  to  fracture  initiation,  load  and 
CMOD  increase  proportionally  in  a  linear  fashion.  The  crack  tip 
is  observed  with  a  traversing  distance  microscope  to  detect  accu¬ 
rately  start  of  fracture  initiation.  Load  to  fracture  initiation 
was  accurately  measured  by  finding  the  points  on  the  load-verus- 
CMOD  curve  at  which  departure  from  linearity  commenced.  After 
fracture  initiation,  load  and  displacement  under  the  load 
adjusted  themselves,  i.e.  increased  or  decreased,  to  keep  the 
CMOD  equal  to  the  CMOD  requested  by  the  program.  At  fracture 
initiation,  the  crack  propagated  in  a  slow  and  stable  manner 
until  it  stopped.  The  location  of  the  new  crack  tip  was  then 
approximately  recorded  using  transparent  graph  paper,  but  later 
accurately  obtained  from  load  versus-CMOD  records  using  compli¬ 
ance  calibration  technique.  This  procedure  was  then  repeated 
until  the  test  was  completed.  Ingraffea  [8],  Ingraffea  and  Arrea 
[9]  used  the  same  testing  procedure  to  test  rock  and  concrete 
beams,  respectively.  Bocca,  Carpinteri,  and  Valenti  [15],  and 
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Ballatore  and  Carpinteri  [16]  later  used  a  similar  experimental 
procedure . 

3 . 7  Typical  Results 

3.7.1  Load,  Displacement,  and  CMOD  Records.  A  typical 
mixed-mode  specimen  (a=l  in,d=6  in),  as  shown  in  Figure  3.1,  was 
tested  under  CMOD  control.  Its  load-versus-CMOD  and  load-versus- 
load-line  displacement  curves  are  shown  in  Figures  3.4  and  3.5. 

The  load-versus-CMOD  curve  follows  a  straight  line  until  the 
load  reaches  a  critical  value.  As  previously  mentioned,  this 
value  corresponds  to  fracture  initiation.  The  load-versus- 
displacement  curve  shows,  however,  a  slightly  nonlinear  relation 
at  small  loads  caused  by  local  nonlinear  deformation  at  load  and 
reaction  points.  Upon  further  increase  of  load,  the  curve 
becomes  linear  up  to  fracture  initiation.  However,  the  predicted 
slope  is  always  greater  than  the  observed  because  of  (i)  nonlin¬ 
ear  effects  of  supports  and  at  loading  point?  and  (ii)  deforma¬ 
tion  in  the  loading  train. 

3.7.2  Crack  Trajectory.  Figure  3.6  shows  an  observed  crack 
trajectory  for  a  mixed-mode  specimen.  Following  fracture  initia¬ 
tion,  the  crack  propagates  along  a  curve  that  changes  curvature 
in  the  upper  part  of  the  beam.  It  approaches  a  vertical  direc¬ 
tion  toward  the  top  of  the  beam  and  does  not  cross  the  load 
plane. 

3.7.3  Observation  of  Fracture  Surfaces.  Figure  3.7  shows 
some  characteristics  of  the  fracture  surfaces  observed  in  the 
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Figure  3.6  Observed  crack  trajectory  for  a  typical  mixed-mode 
specimen  (a  *  1  inch,  d  *  6  inches) 
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Initial  crack  length 


Saw  cut 
Natural  crack 

Perpendicular  curve  showing  where 
crack  stopped  then  initiated 


Almost  parallel  lines  indicating  stable 
and  slow  crack  growth 


Rough  fracture  surface 


Figure  3.7  Fracture  surface  observations  (after  Williams,  Radon  and  Turner  [27]) 
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tests.  The  surface  was  rather  rough  with  almost  parallel  mark¬ 
ings  of  lines  in  the  direction  of  propagation,  indicating  stable 
and  slow  crack  growth.  The  same  observation  was  made  by 
Williams,  Radon  and  Turner  [27].  The  figure  also  shows  a  series 
of  markings  of  curves  that  run  roughly  perpendicular  to  the  crack 
trajectory.  These  curves  indicate  where  the  crack  stopped  and 
re-initiated.  These  curves  provide,  therefore,  an  accurate 
method  to  locate  the  crack  tip  at  each  step  of  propagation.  In 
the  future,  this  additional  information  will  be  used  to  verify 
crack  length  indirectly  computed  from  compliance  calibration. 

3.8  Summary 

The  limited  number  of  published  data  required  that  further 
experiments  be  performed  to  validate  analytical  developments  for 
both  isotropic  and  orthotropic  materials.  These  experiments  are 
investigating  sensitivities  of  load  to  fracture  initiation  and 
crack  trajectory  due  to  length,  inclination  and  location  of 
initial  crack,  and  ratio  and  orientation  of  principal  material 
toughnesses.  Their  purpose  is  to  create  a  uniform  probabilistic 
distribution  of  these  stochastic  variables. 

A  mixed-mode  specimen  was  carefully  designed  to  ensure 
applicability  of  linear  elastic  fracture  mechanics,  ease  of 
preparation  and  testing  procedure  and  other  criteria.  CMOD 
control  was  used  to  ensure  stable  and  slow  crack  growth.  This  is 
because  the  interest  was  not  only  the  initial  fracture  initiation 
load  but  the  entire  fracture  trajectory  and  load  history. 
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Experimental  results  to  date  on  isotropic  plexiglas  include 
records  of  load-versus-CMOD  and  load-versus-load-line  displace¬ 
ment,  as  well  as  crack  trajectories  and  observations  of  fracture 
surfaces.  All  of  these  results,  and  other  material  property 
measurements  also  performed,  will  be  used  in  the  following  sec¬ 
tions  to  validate  both  deterministic  and  probabilistic  analytical 
formulations  for  isotropic  materials. 

Preliminary  work  completed  on  anisotropic  materials 
includes: 

(i)  Fabrication  of  orthotropic  plexiglas  using  a  special 
procedure  to  produce  desireable  ratio  of  principal 
toughnesses . 

(ii)  Measurements  of  elastic  properties  and  direct  tensile 
strength . 

(iii)  Measurements  of  fracture  toughness  in  principal 


directions. 
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4.0  DETERMINISTIC  FRACTURE  VALIDATION 

The  probabilistic  fracture  mechanics  code  PROFRANC  developed 
at  Cornell  is  based  on  methods  of  reliability  analysis  and 
theories  of  deterministic  linear  elastic  fracture  mechanics.  The 
code  accounts  for  the  uncertainty  in  loads,  material  properties, 
crack  geometry  and  other  uncertain  parameters.  The  validation  of 
PROFRANC  requires  the  evaluation  of  fracture  mechanics  theories 
used  in  this  code.  This  evaluation  is  performed  as  follows.  An 
accurate  finite  element  model  of  the  model  specimen  was  prepared. 
The  models  were  calibrated  to  experimental  measurements  to  mini¬ 
mize  uncertainties  in  these  parameters.  An  analysis  was  then 
performed,  and  predicted  deterministic  results  were  compared  with 
actual  measured  results.  These  results  include  load,  displace¬ 
ment,  and  crack  mouth  opening  displacement  records,  fracture  ini¬ 
tiation  load  and  angle,  and  crack  trajectories. 

4.1  Deterministic  Fracture  Analysis. 

Four  tested  specimens  with  varying  initial  crack  length 
(a=l,2,3,4  in)  but  the  same  crack  location  and  inclination  were 
selected  for  deterministic  validation  of  the  code  (Figure  3.1). 
Values  of  geometrical  and  material  parameters  have  been  obtained 
from  measurements  of  these  specimens.  Each  analysis  was  based  on 
average  values  of  measured  material  properties  and  the  observed 
values  of  geometry  parameters  include  location,  length  and  incli¬ 
nation  of  initial  crack,  length,  depth  and  thickness  of  the 
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specimen.  The  material  parameters  include  modulus  of  elasticity 
and  fracture  toughness. 

The  modulus  of  elasticity  was  obtained  from  the  initial 
stiffness  of  observed  load  versus  CMOD.  The  average  value  of  the 
modulus  of  elasticity  for  the  four  selected  specimens  was  474,000 
psi.  The  variation  of  observed  modulus  of  elasticity  about  this 
value  was  6%.  The  observed  average  modulus  of  elasticity  is  very 
close  to  the  value  of  450,000  psi  provided  by  the  manufacturer  in 
Table  3.1. 

The  plane  strain  fracture  toughness  was  obtained  from  a 
single  edge  notch  bending  fracture  toughness  test.  The  average 
value  of  fracture  toughness  was  found  to  be  944  lb-in”3/2.  The 
variation  about  this  value  of  measured  fracture  toughness  was 
1.8%. 

PROFRANC  running  in  a  deterministic  mode  was  used  to  predict 
fracture  initiation  loads,  crack  trajectories,  and  load-versus- 
CMOD  and  displacement.  The  output  is  consistent  with  experimen¬ 
tal  results.  The  analysis  is  based  on  the  maximum  circumferen¬ 
tial  stress  theory.  Several  finite  elements  meshes  of  decreasing 
size  have  been  used  to  insure  computation  accuracy. 

4.2  Experimental  Measurements  and  Deterministic  Predictions 

4.2.1  Load  Versus  CMOD  Records .  Figures  4.1  to  4.4  show 
plots  of  observed  and  predicted  load-versus-CMOD  for  all  four 
specimens.  The  agreement  between  experimental  and  predicted 
results  is  remarkable.  The  post-peak  portion  of  both  curves  show 


I 


39 


Q 

O 

2 

O 

(A 

3 

(A 

CA 

w 

® 

® 

£ 

> 

o 

( 

T3 

c 

CO 

O 

CO 

II 

"O 

® 

II 

•O 

O 

w' 

© 

® 

£ 

E. 

o 

a. 

c 

T3 

C 

CM 

CD 

II 

T3 

CO 

a> 

w 

CM 

3 

CA 

at 

CO 

C 

® 

® 

E 

E 

‘Q 

o 

® 

w 

CL 

CA 

k_ 

Q. 

o 

CM 

£ 

3 

o> 

il 

(•puosnoiu.) 
(spunod)  poo~] 


o 


Q 

O 

2 

O 

i 

CO 

© 

©  £ 

*  £ 

T3  -E 

§  ® 
“  II 

©  ■o 

o  © 
=o  ® 

fl)  ^ 

£  u 
a.  c 


co 

S 

3 

o> 

il 


(spunod)  poo~] 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


41 


TT 

^ r 

® 

k_ 

3 

05 

il 


(•punod)  poo~| 


42 


also  a  very  good  agreement.  This  satisfactory  remarkable  per¬ 
formance  of  PROFRANC  was  expected  because  plexiglas  has  a  very 
small  process  zone  compared  to  specimen  dimensions,  so  that 
linear  elastic  fracture  mechanics  analysis  is  valid. 

Figure  4 . 5  shows  a  plot  of  measured  and  predicted  load- 
displacement  curves.  Prediction  curves  are  stiff er  because  (i) 
nonlinear  effects  at  supports  and  at  loading  point  and  (ii) 
deformation  of  the  loading  table  has  not  been  included  in  the 
numerical  analysis. 

Such  a  disagreement  was  not  observed  in  Figure  4.1  to  4.4 
because  load  versus  CMOD  measurement  is  not  affected  by  either 
nonlinear  effects  nor  rigidity  of  the  system. 

4.2.2  Fracture  Initiation  Load.  Table  4.1  shows  measured 
and  predicted  fracture  initiation  loads  for  all  specimens.  Both 
the  maximum  circumferential  stress  theory  and  the  minimum  strain 
energy  density  theory  were  used  in  predicting  fracture  initiation 
loads . 

Agreement  between  analysis  and  experiment  is  reasonable. 
However,  there  are  differences  between  predictions  of  various 
fracture  mechanics  theories.  These  differences  suggest  another 
source  of  uncertainty  that  is  related  to  current  understanding  of 
fracture  phenomena.  For  example,  fracture  initiation  loads  pre¬ 
dicted  by  the  minimum  strain  energy  density  theory  are  higher 
than  those  obtained  by  the  maximum  circumferential  stress  theory. 
Differences  increase  with  the  ratio  of  mode  II  to  mode  I  stress 
intensity  factors  or  crack  length.  A  similar  observation  was 
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Displacement  (in) 

Figure  4.5  Plots  of  measured  and  predicted  load-versus-load  line-displacement 
for  specimen  #1  (a  =  1  inch,  d  =  6  inches) 
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reported  by  Ewing,  Swedlow  and  Williams  [6]  for  plexiglas  panels 
with  an  inclined  initial  edge  crack  in  tension  and  pure  bending. 

4.2.3  Initial  Fracture  Initiation  Angle.  Fracture  initia¬ 
tion  angles  were  measured  using  an  optical  comparator,  with  a 
magnification  of  10  and  an  accuracy  of  +/”  1/4  degree. 

Table  4.2  shows  a  comparison  between  measured  and  predicted 
fracture  initiation  angles.  Both  the  maximum  circumferential 
stress  theory  and  the  minimum  strain  energy  density  theory  were 
employed  to  predict  the  fracture  initiation  angle  of  each 
specimen. 

The  comparison  indicates  a  reasonable  agreement  between 
experiment  and  predictions. 

4.2.4.  Crack  Trajectory.  Figures  4.6  to  4.9  show  plots  of 
observed  and  predicted  crack  trajectories  for  all  four  specimens. 
The  predicted  trajectories  are  in  good  agreement  with  actual 
trajectories.  However,  there  is  a  consistent  bias  of  the  pre¬ 
dicted  trajectories  toward  the  load  point  in  all  predicted  plots. 
This  is  caused  by  differences  between  modeled  and  actual  boundary 
conditions . 

4 . 3  Summary 

Results  show  that  PROFRANC  running  in  a  deterministic  mode 
provides  an  accurate  technique  for  calculating  load-CMOD  curves, 
fracture  initiation  load,  fracture  initial  angle,  and  crack 
trajectory.  Therefore,  proposed  mechanical  model  and  numerical 
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Figure  4.6  Plots  of  observed  and  predicted  crack  trajectories 
for  specimen  #1  (a  *  1  inch,  d  *  6  inches) 
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Figure  4.7  Plots  of  observed  and  predicted  crack  trajectcnes 
for  specimen  #2  (a  ■  2  inches,  d»6  inches) 
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Figure  4.8  Plots  of  observed  and  predicted  crack  trajectories 
for  specimen  #3  (a  *  3  inches,  d  ■  6  inches) 
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Figure  4.9  Plots  of  observed  and  predicted  crack  trajectories 
for  specimen  #4  (a  ■  4  inches,  d  ■  6  inches) 


methods  of  analysis  and  sound.  Their  use  in  probabilistic  frac 
ture  mechanics  analysis  is  warranted. 
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5.0  EFFECTS  OF  HOLES  ON  MIXED-MODE  FRACTURE 

5 . 1  Introduction 

The  presence  of  holes  in  certain  structures  such  as  aircraft 
is  sometimes  a  requirement  set  by  the  designer.  The  effects  of 
holes  on  cracks  present  in  a  structure  can  be  quite  significant 
in  that  they  can  alter  their  trajectories.  Holes  in  a  structure 
can  be  an  additional  source  of  uncertainty. 

This  section  describes  crack  trajectories  obtained  experi¬ 
mentally  and  predicted  by  PROFRANC  (PRObabilistic  FRacture 
ANalysis  Code)  for  a  structure  of  geometry  shown  in  Fig.  5.1. 

The  following  sections  describe  testing  procedures  and  analyses 
of  these  specimens.  Comparisons  are  also  given  between  crack 
trajectories  of  specimens  with  and  without  holes,  and  between 
observed  and  predicted  results  for  those  with  holes.  Finally,  a 
discussion  on  sensitivity  to  initial  conditions  is  presented. 

5.2  Mixed-Mode  Fracture  Testing  and  Analysis  of  Specimens  with 

Holes 

5.2.1  Experimental  Testing.  Seven  plexiglas  specimens  with 
holes  were  tested  to  study  the  effects  of  holes  on  mixed-mode 
crack  trajectories.  The  overall  dimensions  were  held  constant, 
as  were  the  number,  size,  and  location  of  holes.  The  location, 
d,  and  length,  a,  of  the  starter  crack  were  varied  from  specimen 
to  specimen  as  shown  typically  for  specimens  1  to  3  in  Figures 

5.2  through  5.4,  respectively.  The  holes  were  drilled  very 
slowly  to  reduce  the  likelihood  of  damage  and  residual  stresses 
around  them.  The  loading  of  the  specimens  proceeded  in  a  manner 
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Figure  5.2  Finite  element  mesh  of  specimen  #t  with  holes 
(a  -  1  inch  d  «  6  inches) 
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Figure  5.3  Finite  element  mesh  of  specimen  #2  with  holes 
(a  -  2.5  inches,  d  -  6  inches) 
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Figure  5.4  Finite  element  mesh  of  specimen  #3  with  holes 
(a  •  1  £  inches,  d  -  5  inches) 
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identical  to  that  previously  described  for  specimens  without 
holes. 

5.2.2.  Analysis  of  Specimens.  Typical  meshes  used  for 
PROFRANC  predictions  are  shown  in  Figures  5.2  through  5.4.  A  new 
feature,  automatic  propagation,  was  implemented  in  PROFRANC  to 
decrease  user-induced  bias  in  the  simulations.  With  this  fea¬ 
ture,  the  entire  simulation  process  is  controlled  by  PROFRANC: 

1)  An  analysis  of  the  initial  configuration  is  performed. 

2)  SIF's  are  computed. 

3)  Angle  of  propagation  is  computed  using  the  theory. 

4)  The  crack  is  extended  in  this  direction  by  a  constant 
pre-specified  amount,  here  one-quarter  of  an  inch. 

5)  Automatic  remeshing  is  performed  with  local  refinements 
around  the  crack  tip. 

6)  A  new  analysis  is  performed,  and 

7)  Steps  2  through  6  are  repeated  until  the  crack  tip 

" approaches"  a  hole,  or  the  user  terminates  the  simula¬ 
tion  because  the  length  of  the  predicted  crack  exceeds 
that  observed. 

A  current  limitation  in  PROFRANC  is  that  a  crack  cannot 
intersect  a  free  surface,  such  as  a  hole,  automatically.  In  the 
specimens  simulated  here,  the  predicted  crack  was  propagated 
until  it  "approached"  to  within  about  1/10  of  a  hole  radius  from 
the  edge  of  a  hole. 
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5.3  Effect  of  Holes  on  Mixed-Mode  Crack  Trajectories 

5.3.1  Experimental  Observations.  To  study  the  effects  of 
holes  on  mixed-mode  fracture,  the  crack  trajectories  obtained 
from  specimens  with  and  without  holes  were  compared.  Except  for 
the  holes,  the  specimens  had  nominally  identical  geometrical 
parameters.  Figures  5.5  to  5.7  show  comparison  plots  between 
observed  crack  trajectories  of  specimens  1,  2  and  3  with  and 
without  holes. 

It  is  clear  from  these  plots  that  the  effect  of  any  hole  on 
a  crack  trajectory  becomes  significant  when  the  crack  approaches 
the  hole.  This  is  consistent  with  the  Saint-Venant  principle, 
which  indicates  that  a  hole  begins  to  affect  crack  trajectory 
when  the  crack  tip  is  at  a  distance  of  about  one  hole-diameter 
from  the  hole.  Therefore,  the  hole-diameter  can  be  a  critical 
parameter.  Note  also  that  the  fracture  process  may  be  arrested 
by  a  hole  provided  that  it  is  in  a  particular  location  relative 
to  what  would  have  been  the  trajectory  in  the  absence  of  a  hole, 
although  arrest  does  not  always  occur  in  a  straightforward  man¬ 
ner,  as  indicated  in  Figure  5.7. 

If  the  hole  is  placed  in  front  of  the  crack  trajectory  cor¬ 
responding  to  a  specimen  without  holes,  then  the  crack  might  go 
straight  into  the  hole.  This  is  observed  in  the  first  specimen 
(a  =»  1  in,  d  =  6  in),  as  shown  in  Figure  5.5.  If  the  hole,  on 
the  other  hand,  is  placed  slightly  away  from  this  crack  path,  the 
crack  deviates  from  the  trajectory  of  the  specimen  without  holes 
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Figure  5.5  Plots  of  observed  crack  trajectories  for  specimen  #1 
with  and  without  holes  (a  -  1  inch  d  •  6  inches) 
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Figure  5.6  plots  of  observed  crack  trajectories  for  specimen  #2 
with  and  without  holes  (a  »  2.5  inches,  d  •  6  inches) 
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Distance  from  centerline  of  beam  (inches) 

Figure  5.7  Plots  of  observed  crack  trajectories  for  specimen  #3 
with  and  without  holes  (a  •  1 .5  inches,  d  -  5  inches) 
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and  is  attracted  toward  the  hole,  as  shown  in  Figure  5.6  as  the 
crack  approaches  the  middle  hole. 

5.3.2  Comparison  of  Observed  and  Predicted  Trajectories 

A  deterministic  analysis  of  each  of  the  seven  specimens  with 
holes  was  performed  using  PROFRANC.  A  comparison  between  pre¬ 
dicted  and  observed  crack  trajectories  is  shown  in  Figures  5.8  to 
5.14. 

In  general  deterministic  predictions  were  in  very  good 
agreement  with  observations.  PROFRANC  predicted  arrest  in  the 
correct  hole  in  five  of  the  six  instances  in  which  arrest 
occurred.  Deviation  from  observed  trajectory  was,  with  one 
exception,  always  less  than  one-quarter  inch,  or  a  few  percent  of 
the  current  crack  length. 

The  deterministic  prediction  failed  on  Specimen  3,  Fig. 

5.10,  which  had  perhaps  the  most  complex,  and  unexpected  tra¬ 
jectory.  PROFRANC  failed  to  predict  arrest  into  the  loading  side 
of  the  central  hole,  although  influence  of  this  and  the  lower 
hole  on  predicted  trajectory  are  clearly  evident. 

5.4  Sensitivity  to  Initial  Conditions 

Specimens  3  through  7,  Figures  5.10  through  5.14,  comprise  a 
test  of  sensitivity  to  initial  conditions  in  that  they  span  a 
relatively  small  range  of  initial  crack  locations.  As  the  ini¬ 
tial  crack  location  moves  farther  from  ttw  load  line,  three  dis¬ 
tinct  results  are  observed. 

1.  No  arrest,  specimen  4,  Figure  5.11. 


Figure  5.8  Plots  of  observed  and  predicted  crack  trajectories 
for  specimen  #1  (a  -  1  inch  ,  d  -  6  inches) 
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Figure  5.9  Plots  of  observed  and  predicted  crack  trajectories  (a  -  2.5  in. 
d  -  6  in.) 
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Figure  5.10  Plots  of  observed  and  predicted  crack  trajectories  (a  -  1 .5  in. 
d  •  5  in.) 
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1  Plots  of  observed  and  predicted  crack  trajectories  (a  -  1 .5 
d  -  3.5  in.) 
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Figure  5.1 3  Plots  of  observed  and  predicted  crack  trajectories  (a  -  1 
d  -  4.5  in.) 


Distance  from  centerline  of  beam  (inches) 


Figure  5.14  Plots  of  observed  and  predicted  crack  trajectories  (a  -  1 .5 
d  -  4.75  in.) 
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2.  Arrest  in  the  lower  hole,  but  at  markedly  different 
locations.  Specimens  5,  6,  and  7,  Figures  5.12,  5.13,  and  5.14. 

3.  Arrest  on  the  loading  side  of  the  middle  hole,  Specimen 
3,  Figure  5.10. 

Most  interesting  is  the  sensitivity  shown  in  Specimens  3  and  7, 
Figures  5.10  and  5.14.  The  difference  between  these  specimens 
was  one-quarter  inch  in  initial  crack  location.  This  represents 
less  than  2  percent  of  the  span  of  the  beam.  That  is,  an  uncer¬ 
tainty  of  this  small  order  in  initial  crack  location  can  lead  to 
the  large  difference  in  observed  result.  One  could  conjecture 
that  an  initial  crack  located  between  those  in  these  specimens 
could  produce  a  trajectory  which  does  not  arrest,  like  that  pre¬ 
dicted  by  PROFRANC  for  Specimen  3,  Fig.  5.10. 

Finally,  there  appears  to  be  a  corridor  of  stability  accom¬ 
panying  the  observed  trajectories.  As  long  as  the  prediction 
remains  within  this  corridor,  very  good  comparisons  are  obtained. 
If  the  prediction  deviates  from  this  corridor,  divergence  between 
observation  and  prediction  occurs.  It  is  obvious  that  the  pre¬ 
sence  of  holes  narrows  this  corridor,  and  places  much  more  strin¬ 
gent  demands  on  accuracy  in  numerical  modeling.  The  best  evi¬ 
dence  of  this  is  Specimen  3,  where  a  very  small  difference 
between  observation  and  prediction  as  the  early  stages  was  suffi¬ 
cient,  in  the  presence  of  the  strong  field  gradients  generated  by 
the  holes,  to  cause  divergence. 
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6.0  PROBABILISTIC  FRACTURE  ANALYSIS 

Sources  of  uncertainties  in  the  experimental  specimen 
studied  in  this  report  are: 

(a)  geometrical  parameters: 

(i)  initial  crack  length,  a 

(ii)  distance  of  crack  initiation  point  from  the  load 
line,  d. 

(b)  fracture  toughness,  KIC 

Probabilistic  analyses  are  conducted  in  two  phases.  In 
phase  1,  uncertainties  in  initial  crack  length  and  fracture 
toughness  are  examined.  Phase  2  extends  the  studies  of  phase  1 
by  adding  randomness  in  the  distance  d  of  the  crack  from  the 
center  of  the  beam.  Initial  crack  length  has  been  modeled  as 
uniformly  distributed  over  one  to  four  inches,  i.e.  A  _  U(l, 

4) in.  KIC  has  been  measured  experimentally  from  15  samples.  It 
has  an  average  value  of  944  lb-in”3//2  and  range  (896,968)  lb-in” 
3/^2.  KjC  has  been  modeled  as  a  Gaussian  random  variable  with  the 
above  mean  and  same  coefficient  of  variation  as  for  results  in 
Figs.  2.2  and  2.3,  i.e.  7%.  The  distance,  d,  has  been  modeled  as 
a  uniform  random  variable  with  D  _  U(2,  6) in.  In  the  following 
subsections  the  probability  distributions  of  load  to  fracture 
initiation,  PQ,  and  initial  direction  of  extension,  predicted 

by  the  probabilistic  model  developed  in  [3],  and  implemented  in 
PROFRANC,  are  compared  with  experimental  studies  for  both  the 
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phases.  In  phase  1,  the  predicted  results  have  also  been  com¬ 
pared  with  simulation  studies. 

6.1  Fracture  Initiation 

The  determination  of  probabilities  of  PQ  and  dQ  in  PROFRANC 
involves  several  phases.  First,  the  structure  in  Fig.  6.1  is 
discretized  in  finite  elements  involving  singular  elements  at  the 
crack  tip.  Second,  analyses  are  conducted  to  determine  stress 
intensity  factors  {k11,  kj1}  for  several  values  of  crack  length 
a^,  1=1,  2,  ...  n  in  the  range  of  (1,  4) in  under  unit  load  for 
phase  1.  In  case  of  phase  2,  intensity  factors  {k^^,  k21-1}  are 
determined  for  several  values  of  a^  in  the  above  range  and  crack 
location  d  j ,  j  =1,  2  ...  n  in  the  range  (2,  6) in  under  unit 
load.  Third,  response  surfaces  kj  =  k^a),  1  =  1,  2,  for  phase  1 
and  kj_  ■  k^a,  d)  ,  1  =  1,  2  for  phase  2  are  calibrated  to  the 
resulting  stress  intensity  factors  {k^,  k21}  and  {k^^,  k2i:1} 
respectively.  These  surfaces  give  functional  dependence  between 
stress  intensity  factors  and  random  variables  {A}  and  {A,  D) 
characterizing  initial  crack  geometry.  They  are  used  to  calcu¬ 
late  the  probabilities  of  PQ  and  8Q  based  on  the  theories  of 
linear  elastic  fracture  mechanics. 

According  to  the  maximum  circumferential  tensile  stress 
theory  [2],  the  initial  direction,  9Q  of  crack  growth  is  the 
solution  of  the  equation 

Kj  sin  8  +  Kjj  (3cos(0)  -1)  =  9  (1) 


or 
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X  *  Kj/Kjj  =  (1  -  3  cos (6 )  )/sin(0 )  (2) 

Here  6  is  a  random  variable.  The  distribution  of  8Q  is  obtained 
by  using  equation  (2)  and  by  computing  the  probability  that  8Q  is 
less  than  or  equal  to  a  given  8 .  PROFRANC  has  been  employed  to 
perform  these  computations. 

According  to  the  maximum  circumferential  tensile  theory,  the 
function 

KI,eq(KI,  KII}  "  cos(0/2)  [Kj  cos2  8  -  1.5  KIX  sin  (8) ]  (3) 

is  the  mode  I  stress  intensity  factor  along  a  given  direction  8 
in  a  crack  tip  stress  field  where  Kj  and  Kjj  are  stress  intensity 
factors.  Kx  __  is  maximum  when  8  is  the  solution  of  Equation 

x  t  eq 

(1) .  For  the  structure  under  consideration  Kj  and  KXI  are  stress 

intensities  under  unit  load.  Thus  if  P  is  the  applied  load,  then 

the  crack  initiates  when  P  KT  >  KTO  and  the  load  to  fracture 

i , eq  ic 

initiation,  PQ,  is  given  by  PQ  =  Kic/KI  eq'  Here  pQ  ■’•s  random 
due  to  randomness  in  KIC  and  crack  geometry.  The  distribution  of 
PQ  has  been  obtained  by  PRO FRANC  using  first  and  second  order 
reliability  algorithms  in  the  space  of  random  variables  A  and  KIC 
in  phase  1  and  in  the  space  of  A,  D  and  KIC  in  phase  2. 

The  theoretical  distributions  of  dQ  and  PQ  predicted  by 
PROFRANC  have  been  compared  with  the  results  obtained  from  simu¬ 
lation  and  experiments  in  phase  1.  Figures  6.2  and  6.3  show  the 
distribution  of  SQ  compared  with  the  values  obtained  from  experi¬ 
ments  and  simulation  studies  respectively.  The  probabilistic 
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model  successfully  embraces  all  the  simulation  and  experimental 
values  within  the  predicted  range.  The  distribution  of  8  is 
highly  skewed  due  to  the  nonlinearity  in  the  relation  between 
initiation  angle  and  crack  length  shown  in  Figure  6.4. 

Figures  6.5  and  6.6  show  the  distribution  of  PQ  compared 
with  results  obtained  from  simulation  and  experiments.  Again, 
good  agreement  between  predicted  range  and  simulation  as  well  as 
experimental  results  is  observed.  The  nonlinear  relation  between 
crack  length  and  fracture  initiation  load  (Figure  6.7)  explains 
the  skewness  of  the  density  curve  of  PQ. 

Figure  6.8  shows  probabilities  of  8Q  predicted  by  PROFRANC 
for  fixed  crack  length  (A  =  3  in)  and  uniformly  distributed  dis¬ 
tance  D  in  range  (2,  6)  in.  The  figure  also  shows  experimental 
values  of  8 Q  for  A  =  3  in  and  D  *  2,  3,  4,  5  and  6  in.  The  prob¬ 
abilistic  model  embraces  most  of  the  experimental  values  within 
the  predicted  range.  Figure  6.9  shows  the  probability  of  8 q  pre¬ 
dicted  by  PROFRANC  both  A  and  D  are  random.  From  Figures  6.8  and 
6.9,  the  uncertainty  in  the  initial  direction  of  crack  extension 
increases  when  both  A  and  D  are  random  variables.  This  is  con¬ 
sistent  with  physical  anticipation. 

Figure  6.10  shows  the  distribution  of  PQ  predicted  by 
PROFRANC  for  D  _  U(2,  6)  in,  Klc  _  N(944,  66)  lb-in"3/2  and  A  =  3 
in  and  experimental  values  of  PQ  for  A  =  3  in  and  D  =  2,  3,  4,  5 
and  6  in.  The  agreement  between  predicted  range  and  experimental 
results  is  satisfactory.  Figure  6.11  shows  the  distribution  of 
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Figure  6.5  Probability  density  curve  of  fracture  initiation  load 
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Figure  6.6  Probability  density  curve  of  fracture  initiation  load 


Prodab  1 1 1 ty  dens i ly 


ProbabI 1 i ly  dens! ly 


86 

PQ  predicted  by  PROFRANC  for  A  „„  U(l,  4)  in,  D  „  U(2,  6)  in  and 
Klc  N (944 ,  66)  lb-in-3/2. 

6.2  Hypothesis  Testing 

The  probability  density  curves  obtained  in  phase  1  by 
PROFRANC  for  load  to  fracture  initiation  and  initial  direction  of 
crack  extension  are  quantitatively  validated  in  this  section  by 
Kolmogorov-smirnov  test  and  Wilcoxon's  test.  These  tests  have 
been  discussed  in  section  2  of  the  report.  The  results  from 
Komogorov-Smirnov  for  load  to  fracture  initiation  and  initial 
direction  of  crack  extension  are  given  in  Tables  (6.1)  and  6.2) 
respectively.  The  corresponding  results  by  Wilcoxon's  test  are 
given  in  Tables  (6.3)  and  (6.4). 

6.3  Crack  Trajectory 

Probabilistic  study  of  trajectories  is  based  on  Monte  Carlo 
simulation  and  response  surface  method  in  phase  l  and  only 
response  surface  method  in  phase  2. 

6.3.1  Monte  Carlo  Simulation 

Twenty  eight  samples  of  initial  crack  lengths  are  generated 
for  crack  lengths  uniformly  distributed  from  1  to  4  in  for  tra¬ 
jectory  studies  in  phase  1.  The  distance  of  all  crack  initiation 
locations  from  the  beam  centerline  was  fixed  at  6  in,  (d  =  6  in 
Figure  6.1).  Using  these  samples,  the  structure  has  been  anal¬ 
yzed  by  PROFRANC  in  automatic  crack  propagation  mode.  Thus 
twenty  eight  samples  of  crack  trajectory  was  obtained. 
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Table  6.1  Results  of  analyses  by  Kolmogorov-Smirnov  Test  for  initial 
direction  of  crack  extension.  Random  variable:  Crack  length, 
A  ~  U(1 ,4)in. 
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Table  6.2  Results  of  analyses  by  Kolmogorov-Smirnov  Test  for  load  to 
fracture  initiation.  Random  variables:  Crack  length, 

A  ~  U(1,4)in,  K,c  -  N(944,  66)psi-in<1/2> 
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Table  6.3  Results  of  analyses  by  Wilcoxon  Test  for  initial  direction  of 

crack  extension.  Random  variable:  Crack  length,  A  ~  U(1 ,4)in. 
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Table  6.4  Results  of  analyses  by  Wilcoxon  Test  for  load  to  fracture 
initiation.  Random  variables:  Crack  length,  A  ~  U(1 ,4)in, 
K|C~N(944,  66)psi-in<1/2> 
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Statistical  analysis  was  performed  using  these  samples.  Figure 
6.12  shows  the  sample  trajectories. 

Figure  6.13  shows  the  mean  trajectory  obtained  by  taking  the 
mean  of  the  y-ordinates  of  the  sample  trajectories  at  different 
vertical  sections.  The  vertical  sections  at  values  of  x  greater 
than  -4.2  inch  intersects  with  fewer  than  28  trajectories.  This 
shifts  the  mean  trajectory  towards  the  right  until  it  coincides 
with  the  trajectory  corresponding  to  the  lowest  initial  crack 
length.  The  figure  also  shows  the  trajectories  one  standard 
deviation  away  from  the  mean.  These  trajectories  merge  towards 
the  mean  trajectory  as  the  vertical  sections  are  taken  at  values 
of  x  higher  than  -4.2  inch. 

Figure  6.14  shows  the  cumulative  density  curves  of  trajec¬ 
tories  along  horizontal  sections  at  different  values  of  y.  The 
ordinates  of  these  curves  at  x  give  the  fraction  of  trajectories 
exceeding  y  at  x.  The  spread  of  these  density  curves  decreases 
as  the  horizontal  section  is  raised.  This  is  consistent  with 
physical  anticipation.  In  the  limit  as  y  approaches  infinity  all 
the  trajectories  should  approach  the  load  point.  Thus  at  infin¬ 
ity,  the  density  function  would  be  a  step  function,  located  just 
to  the  left  of  the  load  point. 

6.3.2  Response  Surface 

Phase  1 

A  response  surface  has  been  generated  for  crack  trajectory 
from  sample  crack  trajectories.  The  surface  is  given  by  the 
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Figure  6.13  Second  moment  characteristics  of  sample  trajectories 
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following  function, 

y (A,  X)  =  14.6  +  1.09A  +  3.33X  +  .163A2  +  15x2  -  .0315A3  + 

. 0404A2X  -  . 0037AX2  (4) 

which  defines  the  crack  trajectory  with  respect  to  x  for  a  given 
initial  crack  length,  A.  Figure  6.15  shows  the  trajectories 
obtained  from  Eq.  (4)  compared  with  experimental  trajectories. 

The  response  surface  equation  has  been  used  to  develop  the 
evolution  of  the  trajectory  density  with  x  (Figure  6.16).  The 
density  curves  are  skewed  to  the  left.  This  is  consistent  with 
the  fact  that  the  trajectories  corresponding  to  lower  crack 
lengths  tend  to  converge  towards  those  corresponding  to  higher 
initial  crack  lengths  (Figure  6.12).  This  is  also  observed  from 
Figure  6.17  where  the  curves  indicate  the  y-ordinates  of  trajec¬ 
tories  at  given  values  of  x  with  increasing  crack  length.  The 
slope  of  these  curves  increases  with  initial  crack  length.  Since 
the  initial  crack  length  is  uniformly  distributed,  the  density 
curve  of  y  becomes  skewed  towards  its  lower  values.  As  x 
increases,  the  curves  of  Figure  6.17  become  flatter  and  the  skew¬ 
ness  of  the  density  curves  of  y  in  Figure  6.16  increases.  We 
also  observe  that  the  ordinates  of  the  density  function  in  Figure 
6.17  increase  with  x.  This  is  due  to  the  fact  that  these  density 
curves  account  for  those  trajectories  only  which  are  within  the 
boundary  of  the  beam  along  a  given  vertical  section.  This 
increase  in  the  ordinates  of  the  density  function  with  x  is 
consistent  with  Figure  6.13  where  the  standard  deviation  of  the 
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Figure  6.1 5  Comparison  of  trajectories  from  experiments  and 
response  surface  predictions 
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trajectories  decreases  with  X  shown  by  the  curves  one  standard 
deviation  away  from  the  mean  trajectory. 

Phase  2 

PROFRANC  has  been  applied  in  an  automatic  crack  propagation 
mode  to  determine  crack  trajectories  for  several  values  of  crack 
length  a^,  i  =  1,  2,  ...  n,  in  range  (1,  4) in  and  crack  location 
dj  j  *  1,2,  ...  nf  in  range  (2,  6) in.  A  response  surface  y  = 

j  / 

y(d,a,x)  can  be  calibrated  to  the  resulting  crack  trajectories. 

It  depends  on  random  variables  {A,  D}  specifying  the  initial 
crack  geometry  and  has  the  from 

y  (d,  a,  X)  =  1.3833  +  4.4569d  +  0.057a  +  8.069X  -  0.239d2 

(5) 

+  0.118a2  +  1.296X2  -  0. 0115da  -  0.043ax  +  0.057xd  + 
0.087X3 

for  given  a  and  d.  The  system  of  coordinates  xy  is  shown  in 
Figure  6.1. 

Figure  6.18  shows  histograms  of  crack  trajectories  Y(D,  A, 
x)  for  A  _  U(l,  4) in,  D  _  U(2,  6)in  and  x  =  -5.5,  -4,  -2,  -l  in. 
These  histograms  have  been  obtained  from  the  response  surface  in 
Equation  5  and  probabilistic  characteristics  of  D  and  A  by 
PROFRANC.  They  are  normalized  so  that  each  of  them  has  unit 
area.  There  is  a  significant  difference  between  the  shape  of  the 
histograms  for  various  values  of  x.  For  example,  the  histogram 
for  x  ■  -5.5  in  is  skewed  to  the  left,  a  behavior  discussed  in 


the  previous  report  [31].  On  the  other  hand,  the  histogram  for  x 
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Figure  6.18  Probability  density  curves  for  y-ordinates  of  crack 
trajectory  at  different  vertical  sections  (response 
surface  method) 
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=  -l.o  in  shows  a  much  smaller  spread  and  is  concentrated  at 
higher  values  of  y.  This  is  consistent  with  experimental  and 
analytical  findings  showing  that  differences  between  crack  tra¬ 
jectories  diminish  as  they  approach  the  point  of  application  of 
load. 


6.4  Example  Problem 

Figure  6.19  shows  a  hypothetical  version  of  the  structure 
studied  in  this  report.  In  this  version  a  vertical  crack 
arrester  AB  has  been  used  in  the  structure  at  x  =  -  3.5  in.  The 
important  question  that  one  needs  to  answer  is  the  following: 

"Suppose  one  knows  the  location  of  potential  crack  initia¬ 
tion  and  the  distribution  of  possible  initial  crack  lengths. 

What  is  the  probability  that  the  crack  arrester  will  ever  get  an 
opportunity  of  serving  its  purpose?" 

If  the  distribution  of  initial  crack  lengths  is  uniform  over 
1  in.  to  4  in.  and  the  crack  initiates  at  6  in.  from  the  beam 
centerline  then  from  Figure  6.14  it  is  obvious  that  the  probabil¬ 
ity  that  the  arrester  will  find  a  trajectory  is  zero.  So  the 
arrester  will  never  serve  its  purpose.  If  the  trajectory  were 
placed  at  x  =  -5.0  in,  it  would  capture  trajectories  due  to  any 
of  the  sample  crack  lengths  and  the  probability  that  the  arrester 
would  arrest  the  crack  would  be  unity.  For  arrester  locations 
between  x  =  -3.5  in  and  x  =  -5.0  in,  a  variable  probability  of 
effectiveness  would  be  predicted.  Note  that  these  results  depend 
upon  the  distribution  of  initial  crack  length. 
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6.F  Sensitivity  of  Crack  Trajectories  to  Boundary  Conditions 
Sensitivity  of  crack  trajectories  to  (a)  location  of  load 
point  and  (b) support  conditions  has  been  studied.  This  study  was 
undertaken  to  evaluate  effects  of  differences  between  experiments 
and  computer  models.  For  example,  the  right  support  of  the  struc¬ 
ture  in  Figure  6.1  was  modeled  for  analysis  as  a  roller.  However 
this  boundary  condition  cannot  be  reproduced  in  the  laboratory 
exactly  due  to  friction  and  local  plastic  deformations.  A  hori¬ 
zontal  spring  attached  to  the  roller  used  in  the  analysis  can  be 
used  to  model  these  effects.  The  point  of  application  of  the  load 
remains  on  a  vertical  line  in  the  experiments  but  is  drifting 
laterally  in  the  computer  model  because  it  is  associated  with  a 
node  of  the  finite  element  mesh.  PRO FRANC  has  been  applied  to 
determine  crack  trajectories  for  A  =  1.0  in,  D  =  6.0  in  and 
several  load  locations  c  in  the  range  (5,  7) in  (Figure  6.20).  A 
response  surface,  y  =  y(l,x),  can  then  be  calibrated  to  these 
trajectories,  where  1  =  d  -  c.  It  has  the  expression 

y  =  44.433  +  0.6431  +  25.195X  -  .0232  + 

+  5 . 03  lx2  +  0.1081X  +  0.358X3  (6) 

and  gives  crack  trajectories  for  any  1  in  range  (-1,  l)in. 

Figure  6.21  shows  probabilities  of  crack  trajectories  at 
different  values  of  x  for  A  =  l.Oin,  D  =  6.0in  and  load  location 
c  modeled  as  a  random  variable  uniformly  distributed  over  range 
(5,  7) in.  The  narrow  spread  of  these  curves  indicate  that  the 
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Figure  6.21  Probability  density  curves  for  y-ordinates  of  crack 
trajectory  at  different  vertical  sections  (response 
surface  method) 
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load  position  does  not  have  a  significant  influence  on  the  evolu¬ 
tion  of  cracks  for  this  structure. 

Figure  6.22  shows  the  sensitivity  of  crack  trajectory  to  the 
stiffness  of  a  horizontal  spring  attached  to  the  roller  support 
of  the  beam  for  A  =  1.0 in  and  D  =  6.0  in.  Results  show  that  the 
support  stiffness  has  a  limited  effect  on  crack  trajectory  during 
the  initial  stage  of  propagation.  However,  as  the  crack 
evolves,  the  displacement  at  the  roller  increases  rapidly  and 
compression  caused  by  the  spring  force  becomes  relevant  and 
changes  the  stress  field  at  the  crack  tip.  This  result  in  an 
increase  of  the  ratio  of  mode  II  to  mode  I  stress  intensity  which 
causes  an  increased  deviation  of  predicted  trajectories  from 
experimental  results. 


Support,  stiffness  ■  1QQQ  IP/ln. 


Distance  froa  centerline  of  beaa,  xlln. ) 

Figure  6.22  Crack  trajectories  due  to  various  stiffness  of  the 
horizontal  spring  at  the  roller  support 
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7.0  BIAS  OF  PREDICTED  CRACK  TRAJECTORIES 

It  has  been  observed  that  crack  trajectories  predicted  by 
PROFRANC  using  the  maximum  circumferential  tensile  stress  theory 
deviate  from  experiments.  This  deviation  is  not  significant,  but 
consistently  predicts  trajectories  to  the  same  side  of  those 
observed.  A  similar  bias  was  reported  by  Williams  and  Ewing  [4] 
for  crack  initiation  angle.  The  results  from  the  previous  sec¬ 
tion  appear  to  indicate  that  this  bias  is  not  due  to  systematic 
differences  between  the  load  line  location  or  the  frictional  com¬ 
ponent  of  the  reaction  in  experiments  and  the  model.  It  may  be 
caused  by  nonlinearities  in  the  process  zone  which  are  inconsis¬ 
tent  with  the  maximum  tensile  stress  analysis  based  on  the 
singular  term  of  stress  field  at  crack  tip  [32].  Inclusion  of 
nonsingular  terms  in  the  stress  field  can  significantly 
improve  predictions  of  the  maximum  tensile  stress  theory  [32]. 
Investigation  into  the  cause  of  the  bias  in  predicted  trajec¬ 
tories  is  continuing. 
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8.0  CONCLUSIONS 

A  two-phase  task  was  undertaken  to  validate  PRO FRANC.  Phase 
one  consisted  in  predicting  deterministically  the  outcome  of  a 
subset  of  a  larger  number  of  experiments  in  which  variability  in 
geometry  and  material  parameters  was  purposely  minimized.  The 
purpose  of  this  phase  was  to  verify  that  PROFRANC  could  predict 
nearly  deterministic  events  accurately.  This  phase  was  shown  to 
be  highly  successful. 

This  verification  was  based  on  experimental  results  which 
had  to  be  obtained  within  this  project  due  to  a  paucity  of  com¬ 
prehensive  mixed-mode  fracture  propagation  data  in  the  open 
literature. 

In  Phase  two  all  currently  available  data  involving  inherent 
uncertainties  in  some  material  and  geometrical  parameters  was 
assembled  in  a  probabilistic  framework  and  subsequently  compared 
to  the  probabilistic  predictions  of  PROFRANC  qualitatively  and 
quantitatively.  These  comparisons  were  shown  to  be  very  suc¬ 
cessful.  The  quantitative  comparison  was  performed  by  hypoth¬ 
esis  testing,  which  is  a  mathematical  rule  deciding  whether  to 
except  or  reject  PROFRANC  predictions  using  experimental  data. 

Results  to  date  show  in  particular  that: 

1.  PROFRANC  can  predict  deterministically  load  to  fracture 
initiation,  initial  direction  of  crack  extension,  com¬ 
plete  trajectory  to  failure,  and  various  load- 
displacement  relationships  in  arbitrary  two-dimensional 


structures . 
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2.  PROFRANC  can  accurately  predict  probabilistic  character¬ 
istics  of  fracture  phenomena  for  systems  including 
uncertainties  in  geometrical  and  material  parameters. 
This  predictive  capability  includes  not  only  means  and 
variances  of  all  relevant  fracture-based  responses,  but 
also  their  full  probabilistic  descriptions.  Two  inde¬ 
pendent  hypothesis  tests  namely,  Kolmogorov-Smirnov  and 
wilcoxon  tests,  validated  these  descriptions. 

3.  It  is  possible  to  create  a  controlled,  scientifically 
designed  database  for  both  deterministic  and  probabilis¬ 
tic  response  of  two-dimensional  mixed-mode  fracture 
propagation.  Such  a  database  currently  cannot  be  found 
in  the  open  literature  and  has  been  critical  in  the 
validation  of  PROFRANC. 

4.  A  previously  unanticipated  source  of  uncertainty  was 
uncovered  during  the  experiments  and  simulations: 
variability  in  supposed  boundary  conditions.  This 
source  of  uncertainty  has  been  shown  to  affect  signi¬ 
ficantly  trajectory  of  curvilinear  fracture,  and  can 
have  a  profound  effect  on  the  reliability  of  fail-safe 
design  of  structures  sensitive  to  such  cracking. 
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APPENDIX  A 

MEASUREMENT  OF  ELASTICITY  MODULUS 

A.  1  Introduction 

An  indirect  technique  was  used  to  obtain  an  accurate  measure 
of  Young's  modulus  of  elasticity.  This  technique  uses  data  from 
the  mixed-mode  fracture  tests  and  is  superior  to  a  standard  tensile 
test.  The  following  section  explains  the  details  of  this 
technique . 

A.  2  Description  of  Technique  to  Measure  Young's  Modulus 

A  series  of  deterministic  analyses  of  four  tested  specimens 
with  varying  initial  crack  length  of  nominal  value  (a  =  l,  2,  3, 
4  in)  but  the  same  crack  location  and  inclination  was  first  per¬ 
formed  (Figure  A.l).  Figures  A. 2  to  A. 5  show  the  finite  element 
models  for  all  four  specimens.  In  each  analysis,  a  value  of  unity 
was  used  for  both  Young's  modulus  and  fracture  toughness.  Further¬ 
more,  each  analysis  was  based  on  the  measured  values  of  geometry 
parameters,  namely  location,  length  and  inclination  of  initial 
crack,  length,  depth  and  thickness  of  the  specimen. 

From  each  analysis  the  initial  stiffness  of  the  beam,  defined 
as  the  ratio  of  applied  load  CMOD  applied  load,  is  computed.  The 
actual  initial  stiffness  is  measured  from  the  observed  load-versus- 
CMOD.  This  corresponds  to  the  slope  of  the  load-versus-CMOD  curve 
prior  to  fracture  initiation.  The  ratio  of  observed  stiffness  to 


Figure  A.2  Initial  Finite  element  mesh  of  specimen  #1 
(a  ■  1  inch  ,  d  ■  6  inches) 


Figure  A.3  Initial  Finite  element  mesh  of  specimen  #2 
(a  -  2  inches,  d  -  6  inches) 


Figure  A.4  Initial  Finite  element  mesh  of  specimen  #3 
(a  -  3  inches,  d  -  6  inches) 


Figure  A.5  Initial  Finite  element  mesh  of  specimen  #4 
(a  -  4  inches,  d  -  6  inches) 
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computed  stiffness  equals  the  modulus  of  elasticity  of  the 
specimen . 

Table  A.l  shows  for  each  specimen  the  computed  and  measured 
stiffness  for  each  specimen,  and  the  corresponding  Young's  modulus. 
The  average  value  of  the  modulus  of  elasticity  for  the  four 
selected  specimens  is  474,000  psi  with  a  variation  of  6%  about  this 
value.  The  observed  average  modulus  of  elasticity  is  very  close 
to  the  value  of  450,000  psi  provided  by  the  manufacturer  in  Table 
A. 2. 


Table  A.l.  Computed  arid  Measured  Stiffness  and  Corresponding 
Young's  Modulus  for  Each  Specimen. 
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Specimen 

# 

Initial  Crack 
Length 
(in) 

Computed 
Stiffness 
(lb/ in) 

Measured 

Stiffness 

(lb/in) 

Young 1 s 
Modulus 
(psi) 

1 

1 

0.398 

- T 

186,960 

469,295 

2 

2 

0.184 

85,653 

459,442 

3 

3 

0.097 

47,696 

491,622 

4 

4 

0.053 

_ 1 

25,363 

- L 

475,998 

Note :  Computed  stiffness  is  based  on  a  value  of  1  for  both 

fracture  toughness  and  Young's  modulus. 
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Property 

ASTM 

Method* 

Units 

Plexiglas 

MC 

Mechanical 

Tansila  Strength  (’/•"  specimen* 

0.2VMinuta) 

D638 

Maximum 

PSI 

10.000 

Elongation.  Maximum 

% 

4  5 

Modulus  of  Elasticity 

psi 

450.000 

Flexural  Strength  (Span  depth  ratio  16. 

O.r/Minute) 

D790 

Maximum 

PSI 

15.000 

Deflection.  Maximum 

inches 

0.6 

Modulus  of  Elasticity 

PSI 

450.000 

Compressive  Strength  (0.05'VMinute) 

0695 

Maximum 

PSi 

16.000 

Modulus  of  Elasticity 

PSI 

430.000 

Impact  Strength 

0256 

Charpy  Unnotched 

ft.  lbs./ 

7.0 

W  x  1 "  section 

Rockwell  Hardness 

0785 

— 

M-90* 

Resistance  to  Stress 

Critical  Crazing  Stress 

ARTC  Mod. 

Isopropyl  Alcohol 

Of 

psi 

1.300 

Toluene 

M1L-R-6997 

_ P» 

1.200 

Table  A. 2  Mechanical  properties  of  plexiglas  provided  by  Rhoms 
manufacturing  company 
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APPENDIX  B 

MEASUREMENT  OF  FRACTURE  TOUGHNESS 

B.l  Fracture  Toughness  Testing 

A  standard  plane  strain  fracture  toughness  test,  based  on 
recommendations  of  ASTM  E399-83  [29],  was  performed  to  obtain  the 
fracture  toughness  of  plexiglas.  The  specimen  is  shown  in  Figure 
B.l.  Loading  was  performed  under  CMOD  control  to  assure  a  slow 
stable  crack  growth. 

Readings  of  load  and  CMOD  were  performed  continuously.  At 
fracture  initiation,  the  crack  propagated  in  a  slow  and  stable  man¬ 
ner  until  it  stopped.  The  location  of  the  new  crack  tip  was  then 
approximately  recorded  using  a  transparent  graph  paper.  A  more 
accurate  crack  tip  location  was  later  obtained  from  load-versus- 
CMOD  records  using  the  compliance  calibration  technique.  This 
procedure  was  then  repeated  until  the  test  was  completed. 

From  a  single  specimen,  various  combinations  of  crack  length 
and  corresponding  fracture  initiation  load  were  then  used  to  calcu¬ 
late  the  material  fracture  toughness.  This  is  equivalent  to  test¬ 
ing  several  similar  specimens  with  increasing  initial  crack  length. 

Two  different  methods  were  used  to  calculate  the  fracture 
toughness  of.  The  next  two  sections  explain  in  detail  each  of 
these  methods,  while  the  last  section  presents  a  comparison  between 


them. 
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B.2  Description  of  First  Method 

The  first  method  used  directly  the  measured  crack  tip  loca¬ 
tions  from  the  transparent  graph  paper.  The  corresponding  fracture 
initiation  loads  were  obtained  by  finding  the  points  on  the  load- 
versus-CMOD  curve  at  which  departure  from  linearity  commenced. 

The  apparent  fracture  toughness  of  plexiglas  was  computed 
using  the  equation  proposed  by  Srawley  [30].  The  equation  is  accu¬ 
rate  within  0.5  percent  over  the  entire  range  of  crack  length: 


(B.l) 


where,  a:  initial  crack  length 

p:  fracture  initiation  load 

B:  specimen  thickness 

W:  specimen  width 

S:  specimen  length 

Figure  B.2  shows  a  plot  of  apparent  fracture  toughness  ver¬ 
sus  measured  crack  length.  A  line  fit  through  the  data  would 
give  a  more  or  less  horizontal  line,  which  indicates  that  frac¬ 
ture  toughness  must  be  the  same  at  all  crack  tip  locations.  This 
means  that  the  apparent  fracture  toughness  computed  at  each  crack 
tip  location  is,  indeed,  the  plane  strain  fracture  toughness. 
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Figure  B.2  Plot  of  apparent  fracture  toughness  versus 
measured  crack  length  (First  Method) 
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B.3  Description  of  Second  Method 

A  compliance  calibration  technique  is  used  to  obtain  more 
accurately  the  crack  tip  locations  from  the  load-versus-CMOD 
records . 

The  compliance  calibration  technique  consists  of  obtaining 
from  analysis  a  plot  of  compliance  versus  crack  length.  The  com¬ 
pliance  is  defined  as  the  inverse  of  stiffness.  A  finite  element 
model  of  the  specimen  was  first  prepared,  as  shown  in  Figure  B.3. 
After  performing  an  analysis,  the  crack  tip  location,  and  the 
corresponding  CMOD  and  fracture-initiation  load  were  computed  and 
stored.  The  crack  was  then  propagated,  and  the  procedure  was 
repeated  until  the  crack  reached  the  boundary  prior  to  complete 
failure.  The  compliance  of  the  structure  at  each  crack  tip  loca¬ 
tion  is  defined  as  the  ratio  of  CMOD  to  fracture  initiation  load. 
Figure  B.4  shows  a  plot  of  computed  compliance  versus  crack 
length. 

Since  continuous  readings  of  load-versus-CMOD  were  measured 
during  the  experiments,  the  compliance  corresponding  to  any  crack 
length  could  be  easily  computed  by  dividing  CMOD  by  corresponding 
load.  Knowing  the  compliance  from  experiments,  the  corresponding 
crack  length  can  be  obtained  from  the  compliance  calibration 
curve,  shown  in  Figure  B.4. 

Using  the  same  fracture  initiation  loads  as  in  the  first 
method  and  their  corresponding  CMOD,  the  compliance  at  each  crack 
tip  location  was  computed.  The  corresponding  crack  length  was 
then  obtained  from  the  compliance  calibration  curve.  The 


Figure  B.3  Finite  element  mesh  of  Single-Edge-Notch 
(Bend)  fracture  toughness  specimen 


Figure  B.4  Plot  of  computed  compliance  versus  crack  length 
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apparent  fracture  toughness  was  computed  using  the  same  equation 
used  in  the  first  method.  A  plot  of  apparent  fracture  toughness 
versus  computed  crack  length  is  shown  in  Figure  B.5. 

B.4  Comparison  of  Both  Methods 

Intuitively,  the  second  method  is  more  accurate  than  the 
first  one,  because  it  does  not  include  experimental  errors  made 
in  crack  length  measurement.  In  fact,  the  maximum  relative  error 
between  computed  and  measured  crack  length  was  about  order  of 
11%.  The  second  method  was,  therefore,  used  to  obtain  the  frac¬ 
ture  toughness  of  plexiglas.  The  average  plane-strain 
fracture  toughness  was  944  psi  J in. ,  with  a  range  from  896  to  968 
psi  J in. 
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Figure  B.5  Plot  of  apparent  fracture  toughness  versus 
computed  crack  length  (Second  Method) 
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APPENDIX  C 

DESCRIPTION  OF  AND  USER  MANUAL  FOR 
SOFTWARE  USED  TO  CONTROL  EXPERIMENTS 

C . 1  Introduction 

The  Mixed-mode  Fracture  Experiment  Control  Software  (MIFECS) 
is  an  interactive  menu-driven  program  developed  to  control  testing 
of  mixed-mode  fracture  specimens.  The  software  was  written  in 
QuickBasic  language,  and  is  run  in  the  Microsoft  QuickBasic 
environment . 

C.2  Choice  of  Programming  Environment 

Microsoft  QuickBasic  is  a  powerful  programming  environment  in 
which  applications  can  be  created  for  IBM  personal  computers  and 
compatibles.  It  was  chosen  for  several  reasons  that  include: 

(i)  Support  for  modular  programming  that  enables  one  to 
build  applications  consisting  of  several  module  files. 

(ii)  An  efficient  editor  that  checks  syntax  and  interprets 
each  line  of  code,  and  allows  one  to  run  the  program  at 
any  point  during  the  development  process. 

(iii)  Powerful  debugging  tools  that  provide  efficient  ways  to 
check  the  logic  of  the  program. 

(iv)  A  structured  language  with  many  capabilities  such  as 
subprograms  and  functions,  a  wide  variety  of  structural 
decisions  and  loops,  records,  and  structures,  etc. 


(v)  A  powerful  and  versatile  graphics  package  that  allows 
one  to  write  applications  with  graphics  capabilities  and 
efficient  system-user  interface. 

C.3  Description  of  Software 

MIFECS  consists  of  fifteen  modules.  Each  module  is  assigned 
a  specific  subtask  in  the  program;  for  instance,  the  "DRAW  MENU" 
module  draws  the  menu  page,  and  the  "SAVE_GEO"  module  saves  the 
specimen  geometry  in  an  output  file.  The  modularity  of  the  program 
allows  one  to  (i)  easily  modify  some  parts  of  the  code  without 
changing  the  remaining  parts;  (ii)  easily  expand  the  present  soft¬ 
ware  capabilities;  (iii)  port  the  code  to  other  system  environ¬ 
ments;  and  (iv)  easily  translate  the  code  to  other  languages. 

MIFECS  allows  one  to  read  three  different  input  channels, 
namely  load,  load-line  displacement  and  CMOD.  The  code  also  can 
send  an  output  signal  through  the  CMOD  channel.  The  main  task  of 
the  code  consists  of  simply  increasing  the  CMOD  value  at  constant 
rate  and  reading  the  corresponding  input  channels.  The  same  pro¬ 
cedure  is  then  repeated  until  the  test  is  completed.  At  every 
step,  the  program  waits  for  the  user's  response  to  increase  the 
CMOD  value.  Figure  C.l  shows  a  flowchart  of  MIFECS,  that  explains 
in  more  detail  the  logic  of  the  program. 

MIFECS  uses  a  DT  2818  data  acquisition  system  manufactured  by 
Data  Translation  Company.  The  program  can  be  easily  modified  to 
work  with  other  data  acquisition  systems.  Only  three  modules  need 
to  be  modified:  "INIT  BOARD"  module  which  initializes  the  data 
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Figure  C.1  Flowchart  of  the  program  MIFECS  (Mixed-mode 
Fracture  Experiment  Control  Software) 


35 


TEST  IS  NOT 
COMPLETED 


INCREASE  CMOD 


WAIT  FOR  USER  RESPONSE 
TO  CONTINUE  THE  TEST 
(Hit  <SPACEBAR>  Key) 


READ  INPUT  CHANNELS 
Load,  Displacement,  CMOD 


STORE  DATA  AND  PLOT  IT  ON 
LOAD-DISPLACEMENT  CURVE 


Figure  C.1  Flowchart  of  the  program  MIFECS  (Mixed-mode 
Fracture  Experiment  Control  Software)  (Continued) 
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acquisition  system,  "INCR_CMOD"  module  which  sends  output  signals 
to  increase  CMOD,  "READ_CHAN"  which  reads  all  the  input  channels 
(load,  load-line  displacement  and  CMOD) . 

C.3  User  Manual 

MIFECS  is  an  interactive  menu-driven  program,  whose  objective 
is  to  control  testing  of  mixed-mode  fracture  specimens. 

To  run  the  software,  the  user  must  type  "MIFECS”  at  the  MS- 
DOS  level.  When  the  program  is  activated,  the  main  page,  shown  in 
figure  C.l  appears  on  the  screen.  This  page  consists  of  three  main 
windows.  The  first  window  contains  information  requested  by  the 
program  about  specimen  geometry.  This  includes  specimen  name, 
length,  width  and  thickness,  initial  crack  location,  length  and 
inclination.  The  second  window  contains  a  plot  of  load  versus 
load-line  displacement,  that  indicates  at  which  stage  the  testing 
is.  The  third  window  is  used  to  display  messages  to  the  user  about 
status  of  the  program. 

Prior  to  testing,  the  user  must  input  the  specimen  geometry. 
The  program  then  makes  an  initial  reading  of  all  three  input  chan¬ 
nels  (load,  load-line  displacement  and  CMOD) ,  and  requests  the  user 
to  increase  the  CMOD  value.  The  user  can  do  that  by  striking  the 
<SPACEBAR>  key  on  the  keyboard.  The  code  then  makes  another  set 
of  readings  of  all  input  channels,  and  stores  the  data  in  an  output 
file.  The  procedure  is  then  repeated  until  the  test  is  completed. 
To  quit  the  program  the  user  must  strike  the  <F10>  key. 
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Alter  quitting  the  program,  the  data  is  collected  in  an  output 
file  with  a  name  specified  by  the  user  and  with  extension  ".DAT". 
The  data  can  then  be  analyzed  and  plotted  using  any  spreadsheet 
software,  such  as  LOTUS  or  GRAPHER. 
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APPENDIX  D 

PHOTOGRAPHS  OF  EXPERIMENTAL  SETUP  AND 
OF  A  MIXED-MODE  FRACTURE  TEST 

This  appendix  contains  actual  photographs  of  the  experimental 
setup  and  of  a  typical  mixed-mode  fracture  test  in  progress. 

Figure  D.l  shows  a  photograph  of  the  experimental  setup,  while 
Figure  D.2  shows  the  schematic  of  the  testing  system.  The  configu¬ 
ration  selected  for  testing  is  shown  in  Figure  D.3.  Figure  D.4 
shows  the  specimen  placed  in  the  load  setup.  Figure  D.5  shows  the 
initial  screen  image  of  the  control  program  MIFECS,  which  gets 
displayed  before  starting  the  test.  Figures  D.6  to  D.ll  show 
photographs  of  specimen  and  of  screen  images  obtained  during  a 
typical  mixed-mode  fracture  test.  In  these  figures,  only  three 
stages  are  shown: 

(i)  Just  after  fracture  initiation 

(ii)  When  crack  is  halfway  propagated 

(iii)  When  crack  is  at  its  final  propagation  steps. 


Figure  D.1  Photograph  of  experimental  setup 


Figure  D.5  Photograph  of  initial  screen  image  of  MIFECS 
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Figure  D.7  Photograph  of  screen  image  of  MIFECS  just 
after  fracture  initiation 


Figure  D.8  Photograph  of  specimen  when  crack  is  half  way  propagated 


Figure  D.9  Photograph  of  screen  image  of  MIFECS  when 
crack  is  half  way  propagated 
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Figure  D.10  Photograph  of  specimen  at  final  stage  of  crack  propagation 
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Figure  D.11  Photograph  of  screen  image  of  MIFECS  at  final 
stage  of  crack  propagation 


